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Constraints on the correlation between QSO luminosity 
and host halo mass from high-redshift quasar clustering 
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ABSTRACT 

Recent measurements of high-redshift QSO clustering from the Sloan Digital Sky 
Survey indicate that QSOs at z ~ 4 have a bias (6) ~ 14. We find that this extremely 
high clustering amplitude, combined with the corresponding space density, constrains 
the dispersion in the L—Mh relation to be less than 50% at 99% confidence for the most 
conservative case of a 100% duty cycle. This upper limit to the intrinsic dispersion 
provides as strong a constraint as current upper limits to the intrinsic dispersion in 
the local Mbh — cr relation and the ratio of bolometric to Eddington luminosity of 
luminous QSOs. 

Key words: dark matter - large-scale structure of universe - quasars: general 



1 INTRODUCTION 

It has recently become accepted that quasar activity and 
black hole growth are an integral part of galaxy evolution, 
however a detailed understanding of what triggers quasar 
activity and how they are fueled still eludes us. The lead- 
ing contender for the identity of luminous, high redshift 
QSOs is that they are black holes fe d by by major merg- 
ers of gas-rich g alaxies llCarlbe rg 1990). Recent incarnations 
of such models (Haiman fc Loebll 19981; Cavaliere fc Vittorinil 



a function of l u minos it y and redshift |Croom et al.l [20051; 



Hcnnawi et al 



20061: IPorciani Magliocchctti & Norbcrg 
20041; IPorciani fc Norberd l200d: iMvers et al l l2007al lbllcf 



20001: iKauffmann fc Haehnefd l2000l : Iwvithe fc Loebl 120021 : 



Hopkins et alB oOG) provide a good description of many ob 



served properties of the QSO population. 

The situation is particularly interesting at high redshift, 
where the population of supermassive blac k holes tha t pow- 
ers the QSOs is growing rapidly (see e.g. iFanl [2006, for a 
recent review). To further understand this important phase 
of black hole and galaxy evolution we would like to build a 
model in which QSO activity is tied to the evolving cos- 
mic web of dark matter halos. The relationship between 
QSOs and dark matter halos, their environments and duty 
cycles, can be constrained via obser vations of their space 
density and large-scale clustering |Cole fc Kaised 1 19891 ; 
lHaiman fc HufeoOll ; [Martini fc Weinberdl200ll )." These con- 
straints become particularly sensitive if the QSOs in- 
habit the rarest, most massive halos for w hich the spa- 
tial clustering depends strongly on ha l o mas s (Kaiser 1984; 
Bardeen et ail 1 19861 ; lEfstathiou et all 1 19881 ; ICole fc Kaiserl 



19891 ). 



At redshifts z < 3 the advent of large optical surveys 
for QSOs has led to firm constraints on the clustering as 



Ida Angela et al.l 120081 ; iPadmanabhan et all 120081 ). With 
the Sloan Digital Sky Survey (SDSS) we are now able 
to measure the c lustering of QSOs well e ven at z > 3 
(IShen et al.l |2007|. se e also earlier work by iKundid Il997l ; 
IStephens et alj 19971 ). Interestingly, the correlation length 
of the QSO population increases rapidly with redshift, 
from r = 16.90 ± 1 73 at z ~ 3 to r = 24.30 ± 2.36 at 
z ~ 4. IShen et all (|2007h demonstrate that they can fit 
the observed clus tering and space den s ity of z ~ 4 QSOs 
with the model of iMartini fc Weinberd l|200ll ) provided the 
z ~ 4 QSOs are relatively long lived (i,Q ~ 160 Myr) and 
inhabit halos more massive than ab out 5 x 10 12 h~ 1 MQ. 
For this calculation IShen et al.l (|2007l ) assume that there is 
a monotonic relationship between instantaneous QSO lumi- 
nosity and halo mass, with no scatter. The actual relation 
between instantaneous QSO luminosity and host halo mass 
is expected to include some scatter. Scatter is expected in 
several of the relationships linking the QSO luminosity and 
the host halo mass: the relationship between the host halo 
mass and galaxy bulge, in the relationship between galaxy 
bulge and black hole mass, in the relationship between 
black hole mass and peak luminosity and in the relationship 
between peak and instantaneous QSO luminosity. 

Here we demo nstrate tha t the very high correlation 
length measured bv lShen et alj (|2007l ). when combined with 
the rapid increase in bias for the most massive halos, 
strongly constrains the scatter between instantaneous QSO 



© 0000 RAS 



2 White et al. 




Figure 1. The large-scale bias, (6), vs. the space density for all 
QSOs at z = 4 above a luminosity threshold assuming a log- 
normal relationship between instantaneous QSO luminosity and 
halo mass with dispersion a = 0.1 (long dashed), 0.5 (dashed), 
1 (solid) and 2 (dotted). The case a = (not shown) is almost 
indistinguishable from a = 0.1. In each case the lines run from 
M t = 10 12 to 10 14 - 5 h~ 1 M Q . 



luminosity and halo mass. The essential idea is that any 
scatter increases the contribution from lower mass (and less 
highly biased) halos, so that a measurement of large clus- 
tering amplitude limits the contribution from lower mass 
objects. This constraint is as strong as direct, observational 
constraints at lower redshift on the amount of scatter in the 
various relationships mentioned above. 

We shall work throughout in the ACDM framework and 
adopt the following cosmological parameters: f2 mat = 0.25, 
Qa = 0.75, h — 0.72 and erg = 0.8. Where appropriate we 
shall comment on how our results depend on these particular 
choices. The nex t section out l ines o ur formalism and applies 
it to the data of IShen et al.l {2007), while ^interprets our 
results and discusses future observations. 



2 FORMALISM AND APPLICATION TO 
HIGH- REDSHIFT QSOS 

2.1 The clustering and abundance of QSOs 

We will focus our attention on the space density and large- 
scale bias of QSOs, as these are the easiest to interpret, 
observationally measurable properties of the population. If 
the mean number of QSOs in a halo of mass Mh is N(Mh) 
then their number density and large-scale bias are 



(6) = 



dM h ^-N(M h ) 



dM » ^TT b h{M h )N(M h ) 
aMh 



(1) 

(2) 



where drih /dMh is the (comoving) number density of halos 
per mass interval and bh(Mh) is the bias associated with 
halos of t hat mass. For dn h /dMh we use the fitting func- 
tion from Shct h fc Tormenl (|l999). The large-scale bias is 



slightly more problematic, as our results depend upon mas- 
sive halos being biased and different fits to bh(Mh) from sim- 
ulatio ns have appeared in the literature. ISheth fc Tormenl 
1 19991 ) derived the bias appropriate to their mass func- 
tion fit using the peak background split. A slightly higher 
bias at fixed mass comes fro m assuming ellipsoidal collapse 
l|Sheth. Mo fc Tormenl Il999l ). The larges t bias for halos of 



a give n mass is given by the bias of the Press fc Schechterl 
mass function, as computed bv lCole fc Kaiserl l|l989). 
which is very similar to the fit of Ijing ( 19981 ) for the 
masses of interest. These two forms give biases 30% larger, 
in the ma ss and redshift range o f interest, t han the low- 
est fit (of ISheth fc Tormenl H999I). The fit of ITinker et all 
il2005h lies between that of lSheth. Mo fc Tormenl (Il999l ) and 
ISheth fc Tormenl (l999). For mass-thresholded samples of 
dar k matter halos in a large N-bo dy simulatiorQ], we found 
the iPress fc Schechterl <| 19741 ) and ISheth fc Tormenl (| 19991 ) 
forms bracketed the halo correlation function so this should 
represent the level of current uncertainty. As it provides 
the most conservative estimate of the maximum allowed 
scatter, and the numeri cal data lay closer to the form of 
IPress fc Schechterl (Il974l ). we adopt this as our fiducial bias. 
However, we caution that the statistics in the simulation 
are poor and we also illustrate the effect of choosing the 
ISheth fc Tormenl (|l999l ) form below. In s everal recent simu- 
lation s at other redshifts and masses thelPress & Sch echteil 
dl974h form appears to be a wo rse fit (IShethfc Tormen 
I l999l; ISheth Mo fc Tormenl Il999l ; ISeliak fc Warrenl |2004| ; 
ITinker et al.ll2005l ), but unfortunately there are few direct 
N-body calibrations of bh(M) at the number densities and 
redshifts of interest to us. One could imagine that for ex- 
tre mely rare peaks , who se collapse is much closer to spheri- 
cal (|Bardeen et al.lll986l ), these newer fits a ctually perform 
worse than the origina l formulation (see e.g. ICohn fc White! 
bOOctlDalal et ai1l2008h . This challenging numerical problem 
deserves further investigation. 

To interpret the large-scale clustering measurements of 
IShen et ail l|2007l ) we then must specify N(M h ). We wish to 
choose as simple a model as possible to illustrate the effect 
of scatter on the (b) —n relation, so we assume that the prob- 
ability that a QSO is seen with instantaneous luminosity L 
is log-normally distributed^ around a central value Lo(Mh) 
with a width ox 



P{L\M h )d\nL = 



2ti(Jl 



exp 



In 2 (L/L ) 



2*1 



dlnL 



.(3) 



If a given halo hosts a QSO, the probability that it is above 
some limiting luminosity, i m i n , is then the integral of this 
expression from L m m to infinity, which can be expressed as 
an error function. If we assume that at most a fraction / on of 
halos host active QSOs at any epoch and that Lo oc in 
the region of the turn-on, the mean number of QSOs above 
imin in halos of mass Mh is 



1 The simulation evolved 1024 3 particles in a 1 h 1 Gpc box using 
the Tr eePM code described in lWhitd j|2002h . See lHeitmann et al.l 
l|2007h for a recent comparison with other N-body codes. 

2 We use natural logarithms throughout so that a can be inter- 
preted as a fractional scatter. The log-base-10 based a would be 
ln(10) ~ 2.3 times smaller. 
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N(M h 



-erfc 



InMt/Mh 
V2a 



(4) 



where Lo(Mt) = L m in, cr = ol/oi and we expect a ~ 1. 
We shall assume Eq. Q in what follows and try to con- 
strain a. Figure [T] shows (&(> L)) vs. n(> L)/f on for a 
range of Mt and a. The larger the dispersion the lower 
the bias at fixed space density, as a larger fraction of the 
QSOs are hosted in lower mass halos. Similarly, a decrease 
in / on must be compensated by a decrease in Mt to hold 
n fixed, again leading to a lower (6). In our assumed cos- 
mology the Universe is tv ~ 1.6 Gyr old at z = 4. If we 
assume Lboi/^odd = 0.25 (e.g. iKollmeier et al.ll2006l ) then 
the e-folding time is ~ 0.1 tu, so / on ~ is not un- 

expected. While the inferred lifetimes are then weakly in 
conflict with the upper range o f observational constraints 
(10 6 - 10 s years; iMartinil [2001 ) . these constraints largely 
stem from lower-redshift QSOs. 

Our model is not the most general one describing how 
QSOs inhabit dark matter halos. In particular, there is no 
reason in principle why halos could not host more than one 
QSO. To constrain the functional form of N(Mh) would re- 
quire modeling both the large- and small-scale clustering, 
and additional assumptions about the statistics of QSOs in 
halos. This is beyond the scope of this paper, so we merely 
note that the mass function is very steeply falling in the 
range of halo masses that are of interest. Our statistics are 
thus dominated by the lowest mass halos in our sample. In 
particular, the number density is approximately linear in 
N{Mt). Simply allowing / on > 1 has no impact on (&), but 
it would raise n and slightly weaken our constraints. How- 
ever to increase n significantly above that for / on = 1 would 
require almost every halo of mass Mt to host more than one 
high luminosity QSO which would lead to an unacceptably 
high close pair fraction. For example, if every halo hosted 2 
QSOs the volume averaged correlation function within the 
halo radius would be £ ~ n^rZ 3 . Given the extreme rarity 
of these halos, this would be in conflict with observations 
ilHennawi et~atll2006l ; iMvers et al.ll2007bl . see e.g. Fig. 17 in 
iHopkins et al.ll2008l ). 

At this point it is also easy to see how a change in cos- 
mological parameters affects (b) — ft. For example, if we raise 
08, halos of a fixed mass correspond to less rare peaks. These 
halos are therefore slightly less biased and significantly more 
numerous. To hold ft fixed we would need to increase Mt- 
To hold the bias fixed requires that we increase Mt yet fur- 
ther and therefore to match both the (b) and n the duty 
cycle must increase. For n ~ 10 _7 /i 3 Mpc -3 , a 10% increase 
in as leads to a 4% decrease in (b) , which is insignificant for 
our purposes. In general, for a power-law power spectrum 
with index n e fi, the fractional change in number density at 
fixed bias is 6/(3 + n e ff) times the fractional change in ag. 
Small changes in the other cosmological parameters have 
even smaller effects. 



2.2 Comparison with observations 

This formalism allow s us to use the ob servations of high-z 
QSO clustering from IShen et al.l (|2007) to constrain a and 
/on in our model. To do so we need to estimate th e large-scale 
bias, (b), from the measurements provided in IShen et al.l 
and this requires making several assumptions. 



15 




10 - 
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Figure 2. The large-scale bias, (6), vs. log-normal dispersion, <r, 
for all QSOs brighter than a given L m in chosen such that n = 
5.6 X 10~ 7 h 3 Mpc~ 3 (lower, dotted lines) or 2 X 1CT 7 h 3 Mpc~ 3 
(upper, dashed lines). The shad ed bands ind icate the ±1 and 2<r 
range of bias measured by IShen et al. as described in the 

text. For each triple of lines the upper line assumes / on of Eq. J4j 
is unity while the lower lines assume / on = 0.3 and 0.1. 



First we mu st convert from the clustering quoted by 
IShen et all (120071 ) to large-scale bias. There are several ways 
to do this, and we have chosen to use the amplitude of their 
fits to r~ 2 power-laws over the range 5 < r < 20 ft _1 Mpc. At 
lower z the bias is relatively constant over this range and the 
slope of the dark matter correlation function is close to —2, 
so converting the fit into a measurement of (b) is straightfor- 
ward. However at z > 3 the QSOs are hosted by increasingly 
rare halos, whose bias is becoming more scale dependent. 
Using the mass and halo catalogs from the N-body simula- 
tion described previously we find that the mass correlation 
function is shallower than r~ 2 below 0(10 ft _1 Mpc) but the 
halo correlation function is still quite close to r~ 2 down to 
5/i _1 Mpc. We therefore convert from ro to (b) by matching 



(5) 



at r = 20 ft Mpc, where the bias is close to constant. (This 
is also the scale w h ere w e compared the simulations to the 
iPress fc Schechterl l|l974f ) fitting form to obtain the bias.) 
For convenience w e use the non-linear power spectrum of 
ISmith et al.l (|2003l ) when computing £dm, but we checked 
that this agrees well with the results of the N-body simu- 
lation on the scales of interest and the correction for non- 
linearity is relatively small. 

The clustering measurements of IShen et alj ([2007) are 
averaged over bins in redshift, over which the bias and mass 
correlations are evolving strongly. Fortuitously, the effects 
approximately cancel, leading to a slow evolution in £ qq . We 
estimate the effects of clustering evolution on our constant- 
time constraints by considering 3 models for the evolution: 
passive evolution (which for massive, rare halos corresponds 
almost to £ 9q ^constant), constant bias and constant halo 
mass. Each model predicts t; qq (r, z) f rom which we can com- 
pute the average value measured bv lShen et al.l (|2007l ) as 
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«(r)) 



/ dz (dN/dz) 2 (H/ x 2 Mr,z) 
Jdz (dN/dz) 2 (H/ X 2 ) 



(6) 



where dN/dz is the redshift distribution, _ff is the Hubble 
parameter at redshift z and \ i s the comoving angular di- 
ameter distance t o redshift z. Using dN/dz from Table 1 of 
IShen et all (120071 ) in the two bins 2.9 < z < 3.5 and z > 3.5 
we find that for the passive evolution model the correlation 
length inferred from ({) is within 1% of that inferred from 
£(z = 3) and £(z = 4) respectively. For the constant halo 
mass case £(z) increases with z, but the inferred ro is still 
within 5% of the constant z value for both samples. For the 
constant bias case, £(z) decreases with z, but the inferred 
ro is within 4% of the constant z value for both samples. 
Since the quoted errors on ro are larger than this, we shall 
interpret the quoted correlation lengths as measurements at 

z — 3 and z = 4 respectively. 

Assuming £ 99 = fm/r) 2 . [Shen et all (|2007h quote ro = 
(16.90 ± 1.73) /i _1 Mpc for their "good" sample with 2.9 < 
z < 3.5. Using the dark matter correlation function for our 
cosmology at z = 3 this corresponds to (b) = 7.9 ± 0.8. 
The quoted space density is ft = 5.6 x 10~ 7 /i 3 Mpc- 3 . If we 
use the 'all' sample both ro and hence (b) are 14% lower. 
We will find that this point is relatively unconstraining for 
either choice. 

At z = 4 the best fit is r = (24.30 ± 2.36) /i _1 Mpc, 
again using the "good" sample, corresponding to (6) = 
14.2 ± 1.4, a nd n = W~ r fc 3 Mp c~ 3 . The other clustering fits 
provided by IShen et al l (120071) are discussed below. Note 
however, as pointed out by Ishen et al.l (|2007l ). this space 
density is actual ly an underestimate o f the true space den- 
sity because the lRichards et all (|2006l ) parameterization of 
the LF underestimates the measured space de nsity from 
z ~ 3 -> 4 (see Figure 20 of lRichards et alj|2006l ). At z ~ 4 
this underestimate is approximately a factor of two and we 
therefore use n — 2 x 10 -7 fo 3 Mpc -3 . The higher space den- 
sity leads to a smaller predicted (6), requiring smaller disper- 
sion a between luminosity and halo mass to fit the observed 
clustering (i.e. this provides a more stringent constraint). 

In Fig. [5]we show (b(> L)) as a function of log-normal 
scatter, a, at fixed ft for / on = 1, 0.3 and 0.1. The 68% and 
95% confidence ranges in (b) (above) are also sh own. As can 
be see n, in order to get the large (6) seen by IShen et al.l 
(2007), the scatter at z ~ 4 has to be less than 80% and 
the fraction of halos containing quasars, / on , must be larger 
than ~ 10%. 

The constraints will change if different fits to the clus- 
tering measurements, different fits for bh, or different num- 
ber densities are used. As t he constraints are coming from 
z = 4, we focus on this case. IShen et all (|2007l ) gives 4 fits to 
power- law correlation functions, for "all" or "good" QSO's, 
with fixed power law (ro/r) 2 or varying power law. The fidu- 
cial calculation above was for the "good" QSOs sample fit to 
( r o/r) 2 , which has a x 2 /dof = 0.32. The other combinations 
give: (sample, power-law index, x 2 , 5b/b)= ("all", 2.00, 0.52, 
-15%), ("good", 2.14, 0.32, +6%), ("all", 2.28, 0.50, -6%). 
The central value in Fig. [2]then just moves up and down by 
the shift in bias (the error bars do change slightly in width). 
Our example in Fig. [5] is one of the best \ 2 cases but also 
gives a high bias (and thus is very constraining). 

As mentioned previously, our fiducial bh was chosen to 
give the least stringent constraint. The functional forms 
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Figure 3. The same range of measured bias for z = 4 as shown 
in Fig. [2] but now for two different analytic fits to the theoret- 
ical bias. The Press-Schechter (top, solid) bias fit is as before, 
the Sheth-Tormen (bottom, dashed) is shown below. Again the 
shaded ban ds indicate the ± 1 and 2a range of bias for z = 4 
measured bv lShen et al 1 l|2007h . as described in the text. For each 
triple of lines the upper line assumes f on of Eq. J4]l is unity while 
the lower lines assume f on = 0.3 and 0.1. If the Sheth-Tormen 
bias is used, there is no overlap within the 2cr range of the clus- 
tering measurement. 



of ISheth fc Tormenl ll 19991 ); ISheth. Mo fc Tormenl (|l999h : 
ISeliak fc Warrenl (|2004l )~and iTinker et all (|20TJ5h predict 
lower bh at fixed Mh and hence tighter constraints on 
a and / on . We illustrate the range in Figure [3] which 
compares the re s ults using our fiducial bias to that of 
ISheth fc Tormenl (|l999h . As we discussed before the fits to 
bh(Mh) from N-body simulati ons, where they a r e avai lable, 
tend to lie between those of ISheth fc Tormenl (|l999l ) and 
ISheth. Mo fc Tormenl (|l999l ) which differ from each other 
by 5 — 10% in the mass and redshift range of interest (see 
also lSeliak fc Warren|[200i ; ITinker et al.ll2005l ). Neither the 
Sheth-Tormen or Sheth-Mo-Tormen biases are able to fit the 
measured bias and number density. 

The most unconstraining estimate would be to use 
the lower amplitude of clustering, e.g. the "all" sample, at 
z = 4. We did not choose this because there are a num- 
ber of systematics which could lower the measured clus- 
tering amplitude and the "all" sample has a worse \ 2 ■ 
But if taken in conjunction with the (most conservative) 
Press-Schechter bias and the highest allowed number den- 
sity, the limit can be weakened to allow a — 1 at 95% con- 
fidence for / on = 0.1. It is also possible that QSOs may in- 
habit a special subclass of halos for which the bias is larger 
than the average for that mass. For instance there are in- 
dications that halo history affects clustering in some cases 
Wechsler et "a3200d:IWetzel et al.ll2007l;lGao fc Whitel2007l : 



ICroton. Gao fc Whitell2007l ; |jing, Suto fc Moll2007l ). For the 
very rare and highly biased halos hosting QSOs at z ~ 4 this 
effect may be large enough to weaken our constraint. Unfor- 
tunately the bias of the relevant objects has not yet been 
measured in simulations. 
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3 DISCUSSION 

It appears that the most luminous, highest redshift QSOs 
have instantaneous luminosities which are well correlated 
with their host halo masses. The already strong constraint 
shown here should be improved by measurements of high- 
redshift QSO clustering with yet larger samples. New ob- 
servations will broadly improve the constraint in two ways: 
through wider area observations to identify a larger number 
of QSOs similar in redshift to the SDSS sample, and through 
deeper observations over the same area. Deeper observations 
will identify slightly fainter QSOs at all redshifts and may 
provide a sufficient sample to measure the clustering ampli- 
tude at even higher redshift than z ~ 4, where we expect 
the bias to be even larger. Going wider in area will decrease 
the fractional error on (6) for high-redshift QSOs. As these 
measurements of high-redshift QSO bias will remain Poisson 
limited, the fractional error on (b) will scale as l/Nq ao for 
a survey of iVq SO QSOs. Going deeper to measure luminous 
QSOs at yet higher redshifts is also expected to be useful, 
as the data indicate strong evolution in the clustering am- 
plitude with redshift. 

Figure [4] illustrates how (6) increases with redshift for 
three values of the scatter between QSO luminosity and 
host halo mass assuming f on — 1. At each redshift (b) 
was calculated t o match the Q S O sp ace density evolution 
parametrized bv lRichards et all (|2006l ). although converted 
to our cosmological parameter^ 3 ]. The steepest relation be- 
tween (b) and redshift is for the case of the smallest scatter. 
If /on = 1 and the scatter is minimal, the relative value 
of (b) for different values of the dispersion is a good mea- 
sure of the fractional uncertainty in (6) required to improve 
on our constraint. These ratios demonstrate that the most 
effective constraint from clustering arises from the lowest 
redshift at which / on = 1 and minimal scatter is a reason- 
able approximation. For example, the ratio of (b)(a = 0.1) 
to (6}(cr = 1.0) increases by less than 10% from z = 4 to 
z — 5, while the Nq BO above the SDSS flux limit drops by 
approximately an order of magnitude. Therefore while the 
clustering amplitude increases, the decline in N vo for a flux- 
limited, high-redshift survey more than offsets this gain. By 
contrast, color selection of just very high-redshift (e.g. z > 5) 
QSOs could produce a substantial improvement for a given 
total number of spectra. 

The immediate prospect for improvement in this con- 
straint is completion of the SDS S observatio n s bey ond the 
4000 square degrees employed bv lShen et al.l (|2007l ). These 
observations will simply increase N qso and the error bars 
should scale as N~ s l. On a somewhat longer timescale, the 
proposed Baryon Oscillation Spectroscopic Survey (BOSS3) 
should approximately double the number of z > 3.5 QSOs, 
although it will mostly achieve this by going deeper. Pro- 
vided /on = 1 and small scatter are still reasonable approxi- 
mations to this slightly fainter population (that will still be 
well above the break in the QSO LF), BOSS' factor of two 
improvement in sample size should decrease the error bars 

3 For simplicity we did not correct the Richards et al. (20(3) LF 
parameterization to match their binned LF here. The space den- 
sity has therefore been slightly underestimated from z = 3 — » 4 
and consequently the predicted bias is slightly higher. 

4 http://cosmology.lbl.gov/BOSS 
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Figure 4. (Top) The large-scale bias, (ft), vs. r edshift of halos 
that m atch the integrated QSO space density from lRichards et al] 

(2006) for /on = 1 and dispersion a = 0.1 (long dashed), 0.5 
(dashed), 1 (solid) as in F igure [T] Our bias estimates derived 
from the JShen et alj|2007h estimates are also show n. (Bottom) 
Norm alized number of QSOs per redshift bin from IShen et al.l 

(2007) . The rapid decrease in SDSS QSOs with redshift indicates 
that purely flux-limited surveys at z > 3 can not place as strong 
a constraint on <r at yet higher redshifts, even though those QSOs 
are more strongly biased. 

by approximately a factor of two and produce a powerful im- 
provement in the constraint. For example, if the measured 
bias stayed fixed and the errorbars decreased by a factor of 
two, the upper limit on the observed scatter would improve 
to less than 0.5 (i.e. 0.2 dex) for / on = 1. There are also up- 
coming photometric surveys that go even wider in area and 
deeper: Pan-STARRS^, the Dark Energy Survey (DES3), 
and the Large Synoptic Survey Telescope (LSSlQ). While 
these surveys do not include a dedicated plan for spectro- 
scopic observations of QSOs, they will provide the necessary 
candidate database of these extremely rare objects for spec- 
troscopy. 

On the theoretical side, the calculations presented here 
could also be improved. Better modeling could allow us to 
use a wider range of scales. As large-scale, high-resolution 
simulations become increasingly feasible, it will become pos- 
sible to precisely measure both the bias of all halos and the 
bias as a function of halo history. The scatter between these 
various relationships can also be measured in self-consistent 
numerical sim u lation s of merger-driven black hole growth. 
Irlopkins et al.l (|2007l ) presented a recent analysis of the ob- 
served and simulated scatter between black hole mass and 
various host (bulge) galaxy properties (er, host galaxy mass, 
effective radius) at lower redshifts than we are considering 
here. They find that the intrinsic scatter is generally lower 
than the observed (upper limit) and as small as ~ 0.2 dex 
for combinations that produce a fundamental plane for black 
holes. However, they also note that the strong correlation of 

5 http:/ /pan-starrs. ifa.hawaii.edu 

6 http:/ /www. darkenergysurvey.org 

7 http://www.lsst.org 
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black hole mass with bulge properties implies that the con- 
nection to halo mass is only indirect and dominated by the 
evolution in the typical gas fractions. Measurement of the 
evolu tion of these quantitie s in full cosmological simulations 
(e.g. |Pi Matteo et alj|2008h could provide measurements of 
the expected scatter at the very high redshifts where QSOs 
are observed to be so strongly biased. 

Regardless of these future prospects, the strength of the 
present constraint on scatter in a monotonic relation be- 
tween QSO luminosity and halo mass is already surprising. 
As noted in the introduction, we expect scatter to be present 
because we expect scatter in at least four relationships that 
combine to determine the instantaneous luminosity of a QSO 
in a halo of a given mass: the relationship between instanta- 
neous luminosity and peak luminosity; the relationship be- 
tween peak luminosity and black hole mass; the relationship 
between black hole mass and bulge properties; and the rela- 
tionship between bulge and halo mass. These expectations 
arise from the observed scatter in these individual relations 
at lower redshifts-intriguingly, in all cases the observed scat- 
ter only sets an upper limit on the intrinsic scatter. 

For the relation between black hole mass and 
bulge velocity dispersion known as the Mbh — Q" re - 
lation (IFerrarese fc Merrittl l200d ; iGebhardt et~"ai1 |200d ). 
iTremaine et all (|2002h " showed that the intrinsic dispersion 
is no larger t han 0.25-0 . 3 dex (or 0.58 — 0.70 in our natural 
log notation! . IFerrarese! (|2002h explored the relation between 
Mbh — Mh by using rotation curve data and argued that 
the dispersion between black hole and halo mass could be 
even less than the Mbh — dispersion, although she notes 
that the uncertainty in the transformation between circu- 
lar velocity and halo m ass is not well characterized. Finally, 
iKollmeier et al.l (|2006l ) measured the distribution in Edding- 
ton ratio, iboi/^Edd, for a sample of z = 0.3 — 4 QSOs and 
find the distribution is well described as log-normal with a 
peak at 0.25 and a dispersion of 0.3 dex. As this observed 
dispersion must account for the uncertainty in the mass es- 
timator and bolometric correction in addition to Lboi/iEdd, 
this measurement also sets an upper limit to the intrinsic 
dispersion in I/boi/^Edd- 

In addition to the small scatter in L — Mh, the data also 
suggest that the duty cycle is approximately unity. This may 
be less surprising, given their extremely high luminosities 
and masses of their supermassive black holes. If QSOs are 
limited to accreting at no more than the Eddington rate, 
the luminosities of z > 6 QSOs imply that their central, 
supermassive black holes already exceede d Mbh ~ W 9 M p, 
when the Universe was less than 1 Gyr old l|Fan et al-llioOlb . 
Subsequent near-infrared spectroscopy provides further evi- 
dence that their central black holes a re indeed this massive 
|Barth et alj|2003l : [jiang et all I2007T I. To create such mas- 
sive black holes at early times requires / on ~ 1, modulo the 
uncertainty in the initial seed mass. Together, these strong 
constraints on the duty cycle and scatter in L — Mh for high- 
redshift QSOs provide important new information on how 
supermassive black holes and their host halos grew at early 
times. 

We thank Yue Shen and Cris Porciani for helpful com- 
ments on an earlier draft of this paper. The simulations 
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